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We construct a new perfect action for free staggered fermions, which is more local than the one obtained from 
the standard block average scheme. This pays off in superior properties after a short ranged truncation. This 
action is "gauged by hand" and tested in Schwinger model simulations by means of a new variant of hybrid MC. 
Using "fat links" for the gauge field, we obtain a tiny "pion" mass down to f3 1.5, and the "eta" mass follows 
very closely the prediction of asymptotic scaling. 



It has been observed that a good implementa- 
tion of classically perfect actions provides an ex- 
cellent level of improvement in the sense that lat- 
tice artifacts are drastically suppressed. However, 
the most successful applications have been lim- 
ited to two dimensions so far, and they involve a 
huge number of couplings [QJ^] . This implies that 
corresponding 4d actions are hardly applicable. 
Therefore it is crucial to achieve good improve- 
ment with rather few extra terms, and also this 
can be studied first in d — 2. However, the known 
classically perfect actions of Refs. are not lo- 
cal enough for a very short ranged truncation to 
be sensible, hence - in view of d = 4 - different 
construction methods should be considered. Here 
we present a successful application of "gauging by 
hand" , i.e. we use the truncated perfect couplings 
of the free fermion and insert gauge couplings by 
hand. This leads to a relatively modest overhead 
in the number of couplings compared to the stan- 
dard action, but to an improvement on the same 
level as classical perfection. We conclude that 
this procedure is promising for d = 4. 

Perfect actions are constructed by iterating 
block variable renormalization group transforma- 
tions (RGTs). For massless staggered fermions, 
the RGT should not mix the (pseudo-)flavors in 
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order to preserve the remnant chiral symmetry 
U(l) (8 U(l). This can be achieved by blocking 
each flavor separately As an improvement 
over the standard block averag e (BA) scheme |Q, 
we propose "partial decimation" (PD) ||. The 
free propagator Aba is built by averaging source 
and sink in each block, whereas in the PD scheme 
one only averages the source and treats the sink 
by decimation, or vice versa. Both lead to lo- 
cal perfect actions for free staggered fermions, 
and in both cases we optimize the locality by 
tuning the RGT parameters. The optimization 
criterion is that the mapping down to d = 1, 
A _1 (pi, 0, . . . , 0), only couples nearest neighbors. 
For the PD scheme this is the case for a simple 
6 function blocking at m = ||. However, that 
criterion can be fulfilled at arbitrary mass m. In 
d > 1 the PD scheme leads to a superior degree 
of locality, i.e. to a faster exponential decay of the 
couplings. Hence the truncation can be expected 
to be less harmful. 

The truncation itself should also be done with 
some care. We found mixed periodic bound- 
ary conditions to be optimal: we impose anti- 
periodicity over 6 lattice spacings in the direction 
of odd coupling distance, and periodicity over 6 
lattice spacings in the other direction(s). The re- 
sulting couplings in d = 2 are given in the follow- 
ing Table; for d — 4 we refer again to Ref. §. 
There we also discuss massive fermions, including 
the case of non-degenerated flavors. 
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m — 


BA 


PD 


Cl,0 


0.866299 


0.878234 


Cl,2 


0.066850 


0.060883 


C3.0 


0.016043 


0.010425 


C3,2 


-0.008022 


-0.005212 



The better quality of the truncated perfect ac- 
tion arising from the PD scheme is confirmed by 
the dispersion relation and by the thermodynamic 
scaling behavior |3]||. In the latter case, the scal- 
ing region is extended by one order of magnitude 
compared to the standard action. 



In our application to the Schwinger model 

we obtain a 'V and a "77" particle, since we deal 
with two massless flavors. We use the couplings 
of the PD scheme in the above Table, and we 
include the following lattice paths: 



The crucial question is the acceptance rate; it is 
good down to a minimal value of (3. If we include 
the fat link in ( 1 ) , then it works well down to very 
small ft ([3 H 1.5). Still the simplification in (1) 
accelerates the algorithm by a factor of « 4, while 
the acceptance rate is divided by « 2. 

We generate 3 • 10 5 configurations on a 16 x 16 
lattice, and we measure after multi-steps of 100 
configurations 0. The simulation results ||] are 
unqucnched. 

The spectrum is determined from local bilin- 
ears ||: 

M l _~_{x) = (-l) Xl+X2 feVx+(i,o) -4+(i,<rA] 

We measure the correlators C" 1CT2 (pi, x 2 ) = 
^(M^Si^M^ (0,0))^ . 
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The staple weight u> is tuned so that the small- 
est absolute values of fermion matrix eigenvalues 
are minimized. We find to = 0.238 to be optimal, 
and we use the resulting fat link. 

For the pure gauge part, we use an ultralocal 
plaquette action, which is perfect in d = 2 . 

In our action, a fixed site couples to 24 other 
sites, which is an overhead of a factor of 6 com- 
pared to the standard action. In view of QCD, we 
design an algorithm which avoids a proportional 
increase in compute time: 

Variant of Hybrid Monte Carlo 

(1) Possible Molecular Dynamics steps are iden- 
tified using the standard action. 

(2) The acceptance decision is based on the quasi- 
perfect action. 



• The isovector current correlation is ob- 
tained from C_~. It has only connected 
contributions and yields the "pion" . 

• The isoscalar current correlation is obtained 
from C++. It has disconnected contribu- 
tions too and yields the "eta particle" . 



At p 



3 we evaluate the dispersion relation 
2 1 



E (pi) — m 2 + pf from cosh fits for m w and 
The result is shown in Fig. 1. Then we mea- 
sure these masses with the standard action (from 
Pi = 0) and with the quasi-perfect action (from 
Pi = and p\ — 7r/8), see Fig. 2. The continuum 
predictions are — (scaling), m\ = 2/(tt/3) 
(asymptotic scaling). 

Using the quasi-perfect action described above, 
we find a very small m n , i.e. very good scaling, 
down to [3 ^ 1.5. The rj dispersion follows very 
closely the prediction of asymptotic scaling down 
to (3 w 1. These results — in particular m n , 
which tests the level of improvement — are even 
better than those of Ref. ||, which used trun- 
cated perfect Wilson fermions together with a 

3 The number 100 corresponds to the maximal autocorrela- 
tion time for topological observables (for other quantities 
it is Ri 25). 
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Figure 1. "Meson" dispersion relations at (3 = 3. 

classically perfect vertex function (parameterized 
by 123 independent couplings). 

This shows that a truncated perfect free 
fermion, which is suitably "gauged by hand" , 
can represent a highly improved, short-ranged ac- 
tion. This program is applicable and promising 
for QCD || . 
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Figure 2. The masses and ra, measured at 
various values of (3. 
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